Charge and Current in the Quantum Hall Matrix Model 



T. H. Hansson, J. Kailasvuori, A. Karlhede^ 

' Fysikum, Stockholm University, AlbaNova University Center, SE - 106 91 Stockholm, Sweden 

(Dated: February 2, 2008) 

We extend the quantum Hall matrix model to include couplings to external electric and magnetic 
fields. The associated current suflters from matrix ordering ambiguities even at the classical level. 
We calculate the linear response at low momenta - this is unambigously defined. In particular, 
we obtain the correct fractional quantum Hall conductivity, and the expected density modulations 
in response to a weak and slowly varying magnetic field. These results show that the classical 
quantum Hall matrix models describe important aspects of the dynamics of electrons in the lowest 
Landau level. In the quantum theory the ordering ambiguities are more severe; we discuss possible 
strategies, but we have not been able to construct a good density operator, satisfying the pertinent 
lowest Landau level commutator algebra. 

PACS numbers: PACS numbers: 73.43.Cd, 71.10.Pm 

I. INTRODUCTION 

It is well established that the low-energy effective Lagrangian for a quantum Hall (QH) system is a Chern-Simons 
(CS) gauge theory. For the simplest case of the Laughlin states with filling fraction v = l/(2m+ 1) it takes the form, 

where k = 1/v, and is the external electromagnetic potential. This Lagrangian (and its generalizations to other 
fractions, multi-component systems etc.) correctly describes the quantum Hall conductance, and also the ground 
state degeneracies on higher genus surfaces and edge excitations^. One can also naturally incorporate quasiparticle 
currents and edge modes with appropriate quantum numbers. 

There are, however, important aspects of the QH liquids that are not incorporated in the CS description. For the 
Laughlin fractions described by ^ the quantization of the QH conductance corresponds to quantization of fc, which 
does not follow from any general principle but should be regarded as an input.-O For hierarchical states, the parameter 
k is replaced by a matrix K with quantized entries which again are input parameters. Furthermore, it is difficult, 
within the effective CS framework, to account for the successes of the composite fermion model and to single out the 
experimentally prominent Jain fractions v ~ n/{2pn ± 1)^. 

There have been several recent attempts to formulate effective theories for the composite fcrmions, based on various 
CS mean field theories^ and in these approaches the constraint that keeps the electrons in the lowest Landau level 
(LLL) is very important. On the level of the effective theories this is manifested in the commutation relation 

[pg, pp] = 2z sin (^^{k X p)^ e^^'-'^pj^,, £ = {h/eBf/' , (2) 

for the Fourier components of the LLL-projected density operator and in a dipole form for the quasiparticle 
current operator. Neither of these properties is captured by effective CS theories like (^). For instance, from we 
get p{x) — = ■^^^■'diaj ~ ^b, hence the commutator (0) vanishes [pj^, p^] — —ATTihv{p x k)6{k +p) =0. 

In view of the above, the recently proposed matrix model for the QH effect is intriguing^. This model - which 
can be shown to be a noncommutative version of the effective CS theory - has several appealing features. The 
quantization of k follows from basic principles^ and the specific values corresponding to the Laughlin fractions 
follow by requiring the underlying particles to be fermions, much as in the microscopic Chern-Simons Ginzburg- 
Landau approach^. The classical ground state of the matrix model is an incompressible homogeneous state at the 
Laughlin filling fractions, and the model also reproduces several important aspects of fractional QH physics, such as 
the quasihole charge and statistics, and the presence of edge states. 

In spite of all this the matrix models have so far not added anything to our understanding of QH physics - the 
Laughlin states are very well understood - primarily through the Laughlin wave functions, but also from various 
Chern-Simons mean field theories. The matrix models will be of real physical interest only if they can be generalized 
to describe other states. As discussed above, an important challenge is to understand the Jain fractions, for which 
there are very good wave functions, and where the mean field theory explicitly involves "effective fermions" filling 
auxiliary higher Landau levels. 
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We think that a generahzed matrix model might be useful, since in the existing matrix model it is enough to specify 
that the original particles are fermions in order to get the correct Laughlin fractions. In Laughlin's theory it is the 
combination of fermi statistics and LLL dynamics that determines the wave functions, so, by analogy, we hope that 
the matrix model encodes the correct LLL dynamics. There are already some indications that this is the case. For 
a vanishing non-commutativity parameter 9 it is obviously true by construction, and it has also been shown that a 
certain class of solutions to the finite-dimensional matrix model describes independent point particles moving in a 
stong magnetic field at distances larger than O/l where I is the magnetic lengthX. It is, however, far from clear that 
the matrix model captures the full LLL dynamics in general, and in particular for the incompressible states describing 
the Laughlin liquids. It is the purpose of this paper to study this question, and in particular to construct a conserved 
electromagnetic current that will allow us to calculate response functions and to see if the charge commutator algebra 
(PI is satisfied. 

We proceed by constructing, in the classical matrix model, a particle density and a corresponding conserved current 
density which are then coupled minimally to the external electromagnetic potential. The densities are complicated 
matrix functions, and furthermore, they are not unique - the noncommutativity of the matrices leads to an ordering 
ambiguity already at the classical level. 

Using the infinite dimensional matrix model, which describes the Laughlin states, we calculate the response to weak 
perturbing fields using linear response. This amounts to calculating the ground state density, the QH response and 
the density response to weak and slowly varying magnetic fields. We show that the correct QH results are obtained 
for any choice of matrix ordering in the definition of the densities. 

We then turn to the quantum theory and discuss the construction of the density and current operators. There are 
several constraints on the correct construction: the density must satisfy the commutation relations Q, be properly 
normalized (i.e. reproduce the correct ground state density for homogeneus states), and in the commutative limit 
it must reproduce the physics of non-interacting particles in the LLL. The problem of constructing a density and 
a corresponding conserved current density now involves both a matrix and a quantum ordering problem, and is, in 
general, very involved. 

Our aim is to make this paper fairly self-contained, and acccsiblc to the condensed matter comunity, so we start with 
a short review of the QH matrix models in the next section. The construction of the classical conserved current and 
the corresponding Lagrangian is in section 3, and the calculation of the linear response in section 4. The quantization 
and the discussion of the charge commutation relations is in section 5. We conclude in section 6 with a summary of 
our results, a guess, and some open problems. Some of the material in this paper was presented in^. 

II. THE QUANTUM HALL MATRIX MODEL 

The original version of the matrix model, proposed by Susskind, describing particles of charge — e in a constant 
magnetic field Bz = is given by the Lagrangian^ 

Lo = ^Tr{(x- - aoU) ^abX'' + 20ao} , (3) 

where X°'(t), a = 1,2, and ao(i) are Hermitian matrices ~ the latter is a Lagrange multiplier imposing a constraint. 
The area 9 that enters Lq is the non-commutativity parameter, see (0J below, and in the classical model the average 
density is given by 2Tr9. Quantization introduces h and the magnetic length £ = (h/eBY^^. Since 2tt£^ is the area 
per state in the lowest Landau level, the filling fraction will be related to the dimensionless parameter 9/i^. For 
quantum Hall states, where v is some odd denominator fraction, 9 will be 0{h) and sensitive to quantum corrections. 
In fact, as discussed briefly below, the relation is i^^^ = 9/£'^ + 1. 

The Lagrangian Q is formally very similar to the Lagrangian for one particle moving in a strong perpendicular 
magnetic field Lip = ^eabX°'x^ , where x°'{t) is the position of the particle. In fact, when ^? = 0, © reduces, upon 
solving the constraint, to N copies of Lip {N being the dimension of the matrices X"). We note that, when including 
a potential V{x), the equations of motion obtained from Lip give the drift current i° = -^e'^''dbV{x), but not the 
cyclotron motion. This is as expected since Lip corresponds to infinite cyclotron frequency and hence vanishing 
cyclotron radius for finite velocity. 

As mentioned in the introduction, there is a close connection between the matrix model ^ and the CS theory |^ in 
that the matrix model is equivalent to a non-commutative generalization of the latter. This is how Susskind arrived at 
He noted that in a hydrodynamical description of the QH fiuid one obtains an abelian CS theory with additional 
non-linear terms. He then proposed the non-commutative abelian CS theory - which to lowest nontrivial order in the 
non-commutativity parameter 9 agrees with the non-linear abelian CS theory - as the effective QH theory. 

Varying the Lagrangian (jS)) with respect to the multiplier So yields the constraint 



[X\X^U^i9, 



(4) 
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where the subscript "m" denotes a matrix, as opposed to a quantum, commutator. Taking the trace of Q we see 
that Tr(X^X^) ^ Tr(X^X^) if ^? 7^ 0. Thus the constraint can only be solved by infinite matrices and, as is well known 
from quantum mechanics, the solution is essentially unique. The theory is hence rather trivial - there is only one 
state, corresponding to an incompressible fluid with constant density p — 1/21^9. In order to find other solutions, 
one must modify the constraint Q). This can be done either by introducing external sources^, or by introducing a 
new dynamical field that couples to d^. The latter approach has the advantage of using finite matrices, which means 
that all the ordinary rules for matrix manipulations, such as the cyclicity of the trace, can be used. The finite matrix 
model, due to Polychronakos, is obtained by adding the following piece to the Lagrangian Q 

Lt, ^ ^\ido - a^)^ , (5) 

where 5* is a complex bosonic A^-vector. The X°-^s in Q are now Hermitian N x N matrices and the constraint Q 
is modified into 

[X\x'U=t0-^^^K (6) 
en 

Taking the trace gives 5'^^' = NeB6, which allows for finite TV solutions also when 6 0. 
The Lagrangian L = Lq + Lf, is invariant under the U{N) gauge transformations 

* ^ [/* 

ao UaoU'< - UidoU'' , (7) 

where C/ is a U{N)-matiix. 

This classical model was solved ini. Identifying the eigenvalues of the matrices X^ and X^ as the positions and 
momenta of the N particles respectively, one finds that when the particles are far apart, compared to 9/£, they move 
as N independent particles in a strong perpendicular magnetic field {i.e. governed by Lip = ^eabi°'x^)- On the 
other hand, when a potential cx (X")^ is added - this attracts the particle to the origin - the ground state is indeed 
a circular droplet with constant bulk density p ~ 1/2tt9. The excitation spectrum is that of the Calogero model - in 
particular this means that the low lying excitations can be interpreted as surface modes. 

The most interesting properties of the matrix models are appearant only after quantization. As usual for gauge 
theories, there are two ways to handle the constraint. In the first, which is similar to the Gupta-Bleuler quantization 
of gauge theories, one quantizes the full extended phase space to get the simple canonical commutation relations, 

[-^mm-^rs] "^^^rns^nr (8) 
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and then implement the constraint or ® as a projection operator to define physical states. We shall refer to this 
method as unconstrained quantization. 

Alternatively, one proceeds like in Coloumb gauge quantization of a gauge theory, and first solves the "Gauss' law" 
constraint or ((BJ in some suitably choosen gauge [e.g. X^ diagonal). Writing the Lagrangian in terms of the 
physical variables one can then read off the canonical commutation relations. In this approach, there is no projection 
- all states are physical - but the commutators between the elements of the matrices become complicated. We shall 
refer to this method as constrained quantization. 

The finite matrix model also gives an unambigous relation between the filling fraction v and the Chern-Simons 
level 6 1 i'^ . This was originally derived by Polychronakos by mapping the matrix model onto the quantum Calogero 
model with a coupling constant v^^{v^^ — 1), where v"^ — O/i"^ + I. From the known ground state of the Calogero 
model he could then identify v with the filling fraction^. The quantum shift in the relation between and the 
level O/l"^ is due to the zero point motion of the particles, as explained by HcUcrman and SusskincU^. Assuming the 
original particles to be fermions, one can furthermore show that the filling fraction is quantized to the Laughlin values 
v-'^ = 2m + I. 

We should also mention that invariance under large gauge transformations requires the Chern-Simons level k to be 
an integer independent of the physical interpretation of the theorjiS. For further discussion of the QH matrix models 
we refer to^iLJ^J^. 



III. A GENERALIZED MATRIX MODEL 



We now generalize the matrix model Q and (|SJ| to include coupling to an external electromagnetic field, A^. 
We use the more general finite dimensional model, L — Lq + Li,, which has the additional advantage of allowing 
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cyclic permutations in the trace. At the end of the section we comment on the infinite dimensional case - making 
explicit the modifications this introduces. Since a large constant magnetic field B, perpendicular to the plane, is 
already incorporated in the model, we let 6B denote the magnetic field corresponding to A^^ and assume that it 
has no constant component. To preserve gauge invariance, we must construct a conserved current from the matrix 
variables. It is useful to remember how to treat a single particle moving along the trajectory x{t). Here the particle 
and current densities are given by p{y,t) = 5{y — x{t)) and j{y,t) = x5{y — x{t)), respectively. Current conservation 
V • J ' + dtp — then follows immediately using the chain rule. Using a ^-function to define the particle and current 
densities corresponds to a point particle - but in a non-relativistic theory, one can in fact use any profile function 
fiy — x{t)), corresponding to a rigid charge distribution moving with velocity x{t). 

In the matrix model, the diagonal elements of the matrices, X^^^^, can, in a suitable gauge, be interpreted as the 
guiding center coordinates for the particles when these are far apart compared to 0/^. With this in mind it is natural 
to define the particle density in the matrix model as 

p{y,t) = Tr[Siy'^-X^m , (9) 

where Sly"" — X°') is a matrix-valued kernel. With the same picture in mind, we make the following guess for the 
current density related to the guiding center motion, 

j{y,t) = TriiJC - i[X,aoU)Siy'' ~ X'^m . (10) 

(In the presence of an inhomogeneous magnetic field, there is an additional, purely solenoidal, contribution to the 
current related to the uncancelled cyclotron motion of neighbouring electrons, which will be discussed below.) 

The densities ^ and (|10() are invariant under the L/(A^)-transformations (UJ - provided only that the kernel 
transforms covariantly: 5 — > U6W . However, since are non-commuting objects there are ordering ambiguities 
already at the classical level, leading to inequivalent classical theories. For a kernel of the form 

Siy"" - - / (02 /(fc^y"^ - 6k^) , (11) 

where f{z,x) is dimensionless and analytic in z, UlUI) does indeed give a conserved current. The proof is easy: 



dtTT[ka{y'' - X-)r = -TiJ^My' - X'')rkaX-My- - X')]"-""-' 

m=0 



= -nTikaX^My' - X^)]"-i = -dyaTTX-[kt{y' - X")]" 

and 

9,aTr([X^ao][fc^,(2/''-X^)]") - nTr{[kaX- , ao]My' ~ X'')]"-^ 

= -nTr([fc,X", My' - X^)]"-i]ao) = , (12) 

where we used the cyclic property of the trace. One kernel of the form (|ll|l is 

Sw{y'' -Xn = l^ e^^«to"-^")g(^?fc2) , (13) 

where g{0) = 1. The special case g{x) = 1 is known as Weyl-ordering. The corresponding densities become, in the 
ao = gauge, 

Pk = Tr (e-'^-"^") .g(^^fc2) (14) 
j| = Tr (^X<^e-'''»^''j g{ek^) . 

Note the formal similarity to the one-particle densities. 

However, more general 0(2)-invariant - and J7(iV)-covariant - kernels of the form 

- = / - ^"), ea&fc"(y' - X"), 9k') (15) 
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are possible. (The kernels are assumed to be Hermitian.) For such kernels, 1)10(1 does not generally give a conserved 
current. However, using ^ and current conservation it is straightforward to define a current as an expansion in the 
matrices X"^, k°-{y°- — X"^) and tabk"^ {y'^ ^ X^) that is conserved. An important example of this type is the anti-ordered 
kernel 

Uz -zrz-zt) = J-0^ (.-^),f (.-zt) ^ (,g) 

where z = +'iy'^, k = +ik'^ and Z = X^ + iX'^. It turns out that the classical limit of the charge density operator 
in the quantized matrix model discussed below is given by this kernel. The explicit expression for the corresponding 
conserved current in the — gauge is, using the cyclic property of the trace, 



j{z,z) =jx + ijy = Tr 



^ ' n=Om=0 



(17) 



Having defined a particle density and a conserved current it is straightforward to couple the matrix model Lq + it 
in ^ and Q to an external electromagnetic field by adding 



Unt = ejd^x [p{x, t)Ao{x, t) - j{x, t) ■ A{x, t)] . (18) 

We can cast the resulting Lagrangian in a more transparent form by noting that the current always can be written 
in the form ((l()|l although in general the kernel defining the current will differ from the one defining the charge. For 
example, the kernel defining the current j corresponding to the anti-ordered charge kernel, can be read directly from 
(|17|l replacing the ordinary time derivative with the corresponding covariant derivative. The resulting Lagrangian is, 

L = eTr{(x'^-^[X^ao]m) (^^EabX'' - Aa^ (19) 

+ B0ao + Ao-^5B} + 'f\ido-ao)'f . 

Here the matrices Afj_{X'^, t), /i = 0, 1, 2, are related to the ordinary gauge potential using the kernel 6(jj,y. 

^(X^^)= J d^xA^{x'',t)\^){x'' -X'^), (20) 

and the label fi on the kernel reminds us that it in general depends on which component of it multiplies. In H19|l we 
have also included a term ~ SB(X"-, t), where 6B is related to the magnetic field dB{x°',t), caused by Af^{x°', t), via S 
as in (|20|l . M is the effective mass of the electrons. This is the Simon-Stern-Halperin "magnetic moment" interaction 
that encodes the variation in the Landau energy in a varying magnetic field^**. It gives an additional solenoidal term 
in the current - due to the non-cancellation of the cyclotron motion currents - that is not included in H10|l . 

Note that the charge and current density operators derived from ifT^ do not have an explicit dependence on VE", but 
that there is an implicit dependence since \1/ enters the constraint and thus determines the form of the solutions for 
the matrices X"^. 

It is clear by inspection that the Lagrangian (|19|l is invariant under the noncommutative gauge transformation |(7J), 
and the symmetry under the usual electromagnetic gauge transformation SA^{x) — d^A(x) is ensured by construction 
since the corresponding current is conserved. For general kernels the gauge variation of the matrices A^ is rather 
complicated, but for the special case of Weyl ordering it takes a simple form very reminiscent of the commutative 

22 

case]** 

In the extended model H19|l . the constraint corresponding to © is modified and depends on Aa'. 

[X^e,faX''-|^]™ = 2^0-^^*t^ (21) 
B eB 

and we note that contrary to ^ this does not determine the commutator between X^ and X^. However, if we pick 
the gauge A2 = 0, then we can rewrite the constraint as 

[X\ X' - lii]^ = [X\ i-FV = - 4?**^ ' (22) 
B eB eB 
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where is the momentum conjugate to as calculated from the Lagrangian (|19|l . Clearly, the above construction 
can be carried out in any linear gauge, but we do not know how to handle a general gauge. Note that the simple 
commutation relation (|22(l is between and = eBX^ — eAi, not between X^ and X^, and this makes it hard to 
evaluate the densities © and H1U|) . However, as long as we are interested only in linear response we can expand the 
expression for the density as will be shown in the next section. That the constraint becomes simple when expressed 
in canonically conjugate variables, is in fact very important since it assures that the arguments given in& audi for the 
relation between the parameter 6 and the statistics of the underlying particles hold true also in the extended theory 
(O. 

We now comment on the case of the infinite dimensional matrix model (jJl and its coupling to an electromagnetic 
field. The above analysis used the cyclicity of the trace - e.g. in 1)12(1 and H17|l - which does not hold for infinite 
matrices. Also, defining the particle density as in the finite dimensional case, i.e. by © with the kernel still given 
by H15|) . we can due to the lack of cyclicity no longer use expressions like (|10|) or IjlTI) for the conserved currents. It 
is however straightforward to construct the conserved currents by carefully keeping track of the matrix ordering. So 
will for instance the current corresponding to the Weyl-ordered kernel ((13|l take the form, 



r(f) = / (02^"'^ Tr 



OO n—1 ^ 

E E -,Hk,xT-"'-^x-{~ikdX''r 



_n— m— 



(23) 



Other kernels give similar expressions. It is now clear that we cannot rewrite the interacton Lagrangian (|18|l of the 
generalized infinite matrix model on the simple form (|19|1 . Instead the coupling to the external vector potentials will 
take the more complicated form. 



(24) 



cf. (|17|) and (|18|l . where DqX ~ X — i[X, ciq] is the covariant time derivative. Varying with respect to do one however 
finds that the constraint still has the form H21|) (with ^' = 0), where is given by (|20|l and the kernel (5(^-) is the 
same as for finite matrices. The point is that although trace manipulations using cyclicity are not allowed on the level 
of the Lagrangian, they are permissible in the equations of motion if we (as customary) assume that the variation Sao 
of the multiplier matrix has compact supporlii^. 

In the simplest case when the constraint takes the form [X^,X'^]m = iO, then - and probably any sensible 
kernel H15|) - can be written as a generalized Weyl kernel (|13|) . Also, if one assumes the Hamiltonian to be a scalar 
potential, then by using the equations of motion to replace Xi in H23|l by an expression in Xj's^ one can see that all 
the considered kernels (5(^) are in the infinite case related by factors of the type g{Ok^) with g{0) = 1. This observation 
will be useful in the next section. 

IV. LINEAR RESPONSE IN THE CLASSICAL MATRIX MODEL 

In this section we calculate the response of the infinite dimensional classical matrix model to electromagnetic 
perturbations in the linear approximation. We remind the reader that this model describes only one state with 
constant density p = l/27r0 - nevertheless we obtain non-trivial results. 

We shall consider three quantities, the quantum Hall response to a constant electric field, the ground state density, 
and the response to a weak, static but slowly modulated magnetic field. The two first examples involve response at 
wave vector fc = and the last at small k, i.e. k^9 <^ 1. From the above discussion of the infinite matrix model 
it follows that the classical ordering ambiguities will be of no relevance for these quantities since they amount to 
correction terms ~ k^9. For simplicity, we use the generalized Weyl- kernel H13|) . but any choice of the form Hll|) would 
do. 

A. The quantum Hall response 

We calculate the current response to an applied constant electric field of the form = —EaX°-. To evaluate the 
current, we first need to determine X°- from the equations of motion, a problem very similar to the QH droplet in a 
quadratic matrix potential considered by Polychronakos''. Varying the Lagrangian with respect to X"-, yields22i 

^Xa ^ g'^fc^^i . (25) 
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Because of oc 1, it follows from expressions like H14(l or from the discussion at the end of section III, that the 
electromagnetic current density is given by 

'^ = -^^£ = -"^/'£' (26) 
and for p — v — veB /h a Fourier transformation yields 

J" = -v^eabE'' = -ane'^^'Eb . (27) 
h 

^Prom refsAi we know that if the particles described by the matrix model are fermions, then v^^ — 2m + 1 and (|27|) 
gives the Laughlin values for the quantized Hall conductance an- 



B. The ground state density 

The response to a constant external electric potential SAq will simply give the density of the system 

Br 

p{x) = — = Tr6{x^-X-). (28) 

In the finite case this is difficult to calculate because of the presence of ^I' in the constraint, but in the infinite case 
the calculation can be conveniently done by using coherent states. Defining eigenstates of the lowering operator 

a\a) = a\a) , where a = -^(X^ + iX"^) , (29) 



and using the commutator [a, a^] — 1 appropriate for the infinite matrix model in the absence of external fields, one 
has 

Tr:0(a,at) : = J ^0{a,a) , (30) 

where : : indicates normal ordering of the operator 0(a,a^). 
For the Weyl-ordered kernel (|13|l {g — 1) we now have, 

where k = ki + ik2. Substituting H31(l in H13(l and H28|l we obtain the constant density p = as expected. 

Notice that the factor does not affect the result because of 6^{k). For the same reasons it follows from the 
discussion at the end of section III that the result is independent of the choice of kernel. 



C. Response to a weak perturbing B field 

It is in general hard to calculate the response to an external magnetic field since, as explained above, (|22(l only 
gives an implicit relation for the matrix A'^{X). It can however be done for a weak, and slowly varying perturbing 
field of the form &B{x) — eBsm{q ■ x). 

In the ^ gauge we have A^ — ^ cos((f • x), and the constraint H22|l (with ^' = 0). The annihilation operator 
then becomes 

^ (Xi+z(X2-lii)), (32) 



B 

and kaX"" — 'fj^(ka + fca^) + ^A^- By expanding in e and 9k^, we obtain 

Pfe = Tre-*'="^° =Tr[e-*\/?(^"+'='^')(l-^ii+O(e2))(l + eO(0fc2))]. (33) 
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To evaluate 

5pk = -Tr[e-'Vf + eO{ee) + 0{e^))] , (34) 

is 

we first write 

A^{X) = ±?-[e''?»-^'' +e-''?»x''] ^ ±^[g^yi(9a+gat) +e-»\/|(<?a+<zat) (35) 
2(72 2(72 

and then calculate the trace with the same method as in subsection B above to get 

5pk = -!^^Tr[e-'^('^''^+'^'^^)[e'^(«"+«°*) +e-^^/7(«°+9°*']] + 6i-iO(^^ (36) 

2^2 

— {5^{q- k) + 5^{q + k)) + 0(e) + Oiek") . 
92 V / 

Taking the Fourier transform, we obtain finally 

5p{x) = ^ [sin(g • x) + 0(e) + 0{ek^)\ , (37) 
which is the density response expected in a quantum Hall system. 

V. CHARGE AND CURRENT IN THE QUANTUM MATRIX MODEL 

As discussed in the introduction, the Fourier components of the density operator projected onto the LLL satisfy 
the non-trivial commutation relation Since the hope is that the QH matrix model captures the physics of the 
LLL, we would expect the correctly defined density operator to satisfy (|2J)- First notice that in the infinite matrix 
model, we can easily calculate the poisson brackets of the charge operators 114|) , and turning these into commutators 
we get to lowest order in 9 and fi-, 

[pj:, pp] = if{k X + 0(M, h^) . (38) 

This is at first hand quite surprising, since we would expect the 9 = result to be that of the commutative Chern- 
Simons theory quoted in the introduction, [p^;, pp] = 0. If we however remember that the infinite matrix model only 
describes a single state, with density po = 1/2tt9, we immediately reproduce this result by substituting into the 
RHS of H38|l . Since the infinite matrix model does not allow any density fluctuations, we concentrate on the finite 
model from now on. 

At this point we face a serious operator ordering problem. Note that although the charge as defined in @ with 
a given kernel (|15|) defines a definite ordering of the products of operator-valued matrix elements when the kernel is 
expanded in a power series, any quantum reordering of the operators would correspond to the same classical matrix 
form, so the question of the correct quantum ordering of (|15|) arises. We therefore have ordering ambiguites on two 
levels, on matrix as well as on quantum level. 

If we quantize after solving the constraint, we know the correct answer in the 9 = limit, since we can then use 
commuting diagonal matrices where the entries are simply the coordinates z„ = Xn + iyn and z„ of the N independent 
particles, which satisfy the LLL commutator = 2£'^S,nn- The anti-ordered operator 

N 

p- = Tre-*^e~*^' !z.° ^ e-^^"e-*-"" , (39) 

n=0 

where Z = + iX^, obeys 101 for 6* = 0. This shows that for 9 = the anti-ordered matrix kernel H16() gives 
the correct quantum ordering. Unfortunately, this is not true for 9^0- series expansion shows that p^° no longer 
satisfies 

Below we review an attempt to construct an operator that satisfies the LLL comutator algebra. The basic ob- 
servation is that for iV = 2, the known result (|39|l can be expressed in terms of quadratic operators which can be 
generalized to the 9^0 case but satisfy a ^-independent algebra. In the first version of this paper we erroneously 
claimed to have an expression in these quadratic operators that satisfied the LLL algebra |2Jl. This conclusion was 
due to a logical error as will be pointed out. 
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A. Constrained quantization 



Here we follow Polychronakosi and use the matrices 



- T 



ni^nin 



" ( 1 ^ran ) : 



(40) 



which solve the constraint in the gauge = V eB9{l, 1, . . . 1). Quantizing, the diagonal matrix elements become 
canonically conjugate operators, satisfying (in complex notation), [z^, z„] = 2£'^Smn- Note that 14()|l shows that when 
the particles are far apart \Az\ ^ 0/£, the matrix model describes N independent particles in the lowest Landau level. 
It is also clear that the ordering problem is very hard since Z is a complicated, ^-dependent, operator-valued matrix. 



It is not difficult to find 



operator that satisfies (0) - in fact pk — e 2 '^'■^e 



j^TrZt 



will do. However, this is 



nothing but the center of mass density operator, and will not define a good local particle density except in the trivial 
case of = 1. We shall demand that in addition to satisfying the quantum commutator algebra, the density operator 
must also reduce to the right-hand side of (|39|l in the limit 9 = 0. As advertised above, we have found a solution to 
this problem for 2x2 matrices which is already quite non-trivial, and we now present that construction. 
First we decompose the complex 2x2 matrix Z = Xi + 1X2 as 

Z = CI + zcj'^ + ida'^ 

Z^ = + zcr^ - i'&cr^ , (41) 
where ct* are the Pauli matrices, = 9/{z + z), and where we have introduced center of mass and relative coordinates. 



c 



{zi + Z2) 
[zi - Z2) 



-TtZ 



(42) 



satisfying [(, = [z, z] = while all other commutators vanish. 

For = we expand the matrix exponentials in 139() using the explicit expresions 1)411) and the properties of the 
Pauli matrices to get 



Pk \e=o 



Tre-^'=^e-5'=^^| 



m,n— 



( —ik\2m ( —ik\2n 

(2m)! {2n)\ 



(z2)™(z2)" 



(z^)"'zz(z^)" 
~4~(2m + l)(2n+l) 



Note that since the particles are identical only even powers of the relative coordinate can appear in (|43|l . 
Next define the quadratic operators. 



A; 



iTrZ^ 



B 



^Tr(ZT 
iTrZZt 



{\TrZY = 
- i(TrZt)2 



iTrZTrZt 

4 



z'~§' 
z^ - if 
zz + d'^ . 



It is easy to verify that they commute with C, and satisfy the following ^-independent algebra, 

[CCl = 

[A,B] = 2fA 

[A,B] = -21^ A 

[A, A] = AfB~2£'^ . 



This allows us to define a ^-dependent density operator by substituting 
expansion (|43ll to get. 



(43) 



(44) 



A, z2 



A and zz 



(45) 

B in the series 



' -ik \ 2m I —ik \2n 



(-ikVZm ( 



= 2 



^ (2m)! 

m,n— 



(2n) 



A"'A'' 



4 (2m-t- l)(2n-f 1) 



3 2 ^ 



(46) 



By construction the commutator algebra of the pk is ^-independent, so the commutator ^ can in a ^-independent 
way be reduced to a sum of terms of the form . We can however not reduce these terms to sums of terms of 
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the type A^BA^ or A"^A'^ without using a ^-dependent relation, and thus not conclude that since H43() satisfies the 
algebra l|2Jl, so does H46II . This was the logical error referred to above. 

We now prove that the classical limit of pk in 146|) is the anti-ordered classical operator in 1)39(1 also for finite 6. A 
general 2x2 matrix Z with commuting elements satisfies: 

= - ZTtZ + i {T?Z - TtZ^) 1 = {Z - (i f - Al (47) 

where C and A are defined in 1421) and (|44|l . Using this relation, we can rewrite the classical anti-ordered operator, 

- Tre-^^e-^^' = Tre"* ^e"* (^-^i)e-^(^'-Ci)e-^f . (48) 

After expanding the matrix exponentials and using (|47|) to rewrite all even powers of {Z — <^1)^{Z'^ — CI)" in terms 
of A, A and B, it is a matter of working out a few traces to recover the expansion H4t)|) . 

B. Unconstrained quantization 

In this approach the matrix elements of satisfy the quantum commutation relations (O, and the parameter 9 
enters only via the U{1) part of the constraint © that defines the physical states. Since the density operator is an 
observable it must be U{N) invariant, and since the commutation relations © preserve this invariance, it follows that 
the commutator of two densities must again be U{N) invariant. (Note that since (|46|l is expressed only in traces it 
is a manifestly U{N) invariant expression, althouth it was derived in a particular gauge.) This would seem to imply 
that there could be no 0-dependence in the operators, but this conclusion does not necessarily follow since we are 
considering operators on an extended phase space. The commutator ^ might have additional terms proportional to 
the constraint ©, and this would still give the correct algebra on the physical subspace. However, since the gauge 
constraint explicitly involves this boundary field 4", the density operators can have a non-trivial ^-dependence only if 
they also depend exlicitly on 5", and this would bring us out of the class of operators defined by quantum reorderings 
of @. Although we shall not consider this possibility in this paper, it cannot be excluded a priori. 

Thus restricting ourselves to U{N) invariant combinations of the matrices Z and Z^, satisfying [Zki, Zmn] — 0, 
[ZkijZj^^] — 2£'^SknSim, we note that the classical derivation of (|46|l based on (|48() holds true also in the quantum 
theory since it did not involve any reordering of Z^i^s and Zj^^'s. (In the constrained quantization scheme this 
derivation brakes down since elements of the matrix Z do not commute with each other.) The operators A, A, B and 
C are here again defined in terms of the matrices by H42|) and (|44|l . but they do not satisfy the algebra H45|l . Remember 
that the constrained and unconstrained Z are not unitarily equivalent - in that case they would have to satisfy the 
same algebra - but related via a projection onto the subspace defined by the constraint. The new algebra is however 
very similar to 1)45(1. the only difference being the commutator 

[A, A] = 4^2^-6^^ (49) 

This difference can be absorbed by a quantum reordering of 1(46(1 in which B is replaced by 

S = jTiiZZ^ + Z^Z) - ^TtZ^TtZ = B . (50) 

This gives [A,A\ — M^B — 2£^, while leaving the other commutators unchanged (except for B B). Since the 
operator B has the form of a trace of the number operator, it is not surprising that it differs in the two schemes since 
the number of degrees of freedom are N and N'^ respectively. Note, however, that if the gauge-fixed Z^s are used, 
then B = B, so l(46() with B works in both quantization schemes. 

C. The conserved quantum current 

Just as in the classical case, the current density can be constructed from the particle density by requiering the 
current to be conserved, 

- -77 jk - -77 jk ^ Pk^-7 [Pk,H] , (51) 
2 2 in 

where H is the full quantum Hamiltonian. From this we can derive two expressions for the current. The first is 
obtained simply by evaluating the time derivative of ((46(1 and noting that just as in the derivation of a classical 
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current like (|23|l . a fc or a fc can always be factored from the corresponding expressions, thus defining the complex 
components of the current. The explicit expressions, containing Z and Z^, are not very illuminating. Alternatively, 
we can assume some particular form for the Hamiltonian, evaluate the commutator, and derive an explicit expression 
for the current in terms of the basic matrix operators. We examplify with a Hamiltonian consisting only of an external 
scalar potental, i.e. 

H = V{p)^ I ^V{k,k)p,. (52) 

With this choice, (|51|l takes the form 

P'^-hlw? = ^ y E ^-i^vipMikpT - [kpTW^. , (53) 

where we used Q for the commutator and expanded the exponential and sine functions in a power series. Comparing 
with H51I) we can now directly read off the current, 



which can be written more compactly in real space. A Fourier transformation gives, 

j(^, ^) = T E d:~\p{z, -z)d,V{z, -z)] . (55) 



h 

n=l 



This expression was derived earlier by Martinez and Stone^^ using a second quantized formalism. The above derivation 
is more general in that it applies to any density operator satisfying the algebra l(2Jl, and in particular to the matrix 
model density operator (gHJ- 



VI. SUMMARY AND OPEN PROBLEMS 



We have extended the classical quantum Hall matrix models of Susskind and Polychronakos to include couplings 
to an external field. We gave a general prescription for the construction of a conserved current and noticed that it 
suffers from ordering ambiguities even at the classical level. Nevertheless, zero and low momentum observables can 
be reliably calculated using the infinite matrix model and in particular we showed that the ground state density, the 
quantum Hall response and the response to a weak and slowly varying magnetic field agree with what is expected for 
a fractional QH system. 

It is much harder to construct the correct charge and current operators in the quantized matrix model. The 
difficulties are due to the complicated intermingling of matrix and quantum ordering. In this connection, we also 
notice the difficulties encountered in a recent attempt by Karabali and Sakita to find good particle coordinates in the 
quantum matrix modeUi. It turns out that the most natural guess, i.e. defining the coordinates via the coherent 
states of the annihilation operator Zij, is not correct in that it does not reproduce the Laughlin wave functions. It is 
not unlikely that this is related to the failure of the anti-ordered matrix expression (|16|l to satisfy the correct density 
operator algebra. 

Although not conclusive, our study supports the conjecture that the classical quantum Hall matrix model captures 
important aspects of the LLL physics not present in the usual effective Chern-Simons theory. At the practical level, the 
conclusion is, however, rather disappointing since local observables like charges and currents are rather complicated 
and thus cumbersome to work with. This might perhaps have been anticipated, since the gauge invariant objects in 
non-commutative theories are inherently nonlocal. On the other hand, it is the very fact that a non-commutative 
field theory can be viewed as a commutative theory with an infinite number of higher derivative terms, that makes it 
at all possible that the Landau level structure could be incorporated at the level of an effective Lagrangian. 

We end with a comment on the relation between the classical and the quantum models. It is striking that the 
classical matrix model incorporates features of the FQHE which are related to the Landau level structure and not 
easily described in conventional mean field approaches. The ^-dependent repulsion between the particles is a phase 
space exclusion effect which can be thought of as the classical counterpart of fermion, or more generally anyon, 
statistics. In this way the classical matrix model can emulate a fermionic system with Haldane type pseudopotentals. 
In the case of a quadratic confining potential, this idea is made manifest via the mapping onto the Calogero modei^, 
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and it would be interesting if a more general connection could be made between the QH matrix model and classical 
exclusion statisticaiS,. The analogy between the classical matrix model and the QH system goes even further in that 
the classical counterpart to the density operator algebra |(5J) can be interpreted as a quantum commutator algebra of 
fractionally charged quasi-holesi^. We hope to return to some of these question in the future. 
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There are, however, microscopic derivations of Q, relying on various mean field approximations, where the quantization of 
k follows from a statistical transmutation—. 

This choice, which we here adhere to, results in the one-particle quantum commutator [x, y] = +i£^ and states (j)„ ~ 
z" exp —^\z\^ . More common in the condensed matter community is to choose Bz = — -B, in which case [x,y] = ~i£^ and 
exp-jiyl^p. 

Using II13II . SA^{x) — 9pA(a;) induces the following gauge variations in the matrices 5A^„ = ~ IdA/dX"'] = 

V / mn 

—dTrA/dX^m and SAo = doA where the matrix-valued function A(X°',t) is related to A(x",t) using the Weyl kernel 
lO in lUnj, and the derivative satisfies de'^^"^" /dX'' = -ikae-^''''^\ 

Here we again the cyclicity of the trace in the equations of motion. This assumption can be avoided for a kernel that behaves 
like a (5— function: J &^xx'^5{x,X) — X" - then we have Aq = ~EaX°'. Substituting this into the Lagrangian and varying 
with respect to gives 12511 . 



